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General relativistic deflection of liglit by mass, dipole, and quadrupole moments of gravitational field 
of a moving massive planet in the Solar system is derived in the approximation of the linearized 
Einstein equations. All terms of order 1 /ias are taken into account, parametrized, and classified 
in accordance with their physical origin. The monopolar light-ray deflection, modulated by the 
radial Doppler effect, is associated with the total mass and radial velocity of the gravitating body. 
It displaces the apparent positions of stars in the sky plane radially away from the origin of the 
celestail coordinates associated with the planet. The dipolar deflection of light is due to a trans- 
lational mismatch of the center of mass of the planet and the origin of the planetary coordinates 
caused by the inaccuracy of planetary ephemeris. It can also originate from the difference between 
the null cone for light and that for gravity that is not allowed in general relativity but can exist in 
some of the alternative theories of gravity. The dipolar gravity field pulls the apparent position of 
a star in the plane of the sky in both radial and orthoradial directions with respect to the origin 
of the coordinates. The quadrupolar deflection of light is caused by the physical oblateness, J2 , of 
the planet, but in any practical experiment it will have an admixture of the translation-dependent 
quadrupole due to inaccuracy of planetary ephemeris. This leads to a bias in the estimated value of 
J2 that should be minimized by applying an iterative data reduction method designed to disentangle 
the different multipole moments and to fit out the translation-dependent dipolar and quadrupolar 
components of light defiection. The method of microarcsecond interferometric astrometry has the 
potential of greatly improving the planetary ephemerides, getting unbiased measurements of plan- 
etary quadrupoles, and of thoroughly testing the null-cone structure of the gravitational field and 
the speed of its propagation in the near-zone of the Solar system. We calculate the instantaneous 
patterns of the light-ray defiections caused by the monopole, the dipole and the quadrupole mo- 
ments, and derive equations describing apparent motion of the deflected position of the star in the 
sky plane as the impact parameter of the light ray with respect to the planet changes due to its 
orbital motion. We discuss the observational capabilities of the near-future optical (SIM) and radio 
(SKA) interferometers for detecting the Doppler modulation of the radial deflection, and the dipolar 
and quadrupolar light-ray bendings by the Jupiter and the Saturn. 
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I. INTRODUCTION 

Attaining the level of a microarcsecond (^as) positional accuracy and better will completely revolutionize funda- 
mental astromctry by merging it with rclativistic gravitational physics. Beyond the microarcsecond threshold, one 
will be able to observe a new range of celestial physical phenomena caused by gravitational waves from the early 
universe and various localized astronomical sources, space-time topological defects, moving gravitational lenses, time 
variability of gravitational fields of super-massive binary black holes located in quasars, and many others [1-3]. Fur- 
thermore, this will allow us to test general theory of relativity in the Solar system in a dynamic regime, that is when 
the velocity- and acceleration-dependent components of gravitational field (the metric tensor) of the Sun and planets 
bring about observable relativistic effects in the light deflection, time delay and frequency [4], to an unparalleled 
degree of precision. 

Preliminary calculations [5] reveal that the major planets of the Solar system arc sufficiently massive to pull photons 
by their gravitational fields, which have significant multipolar structures [6], in contrast with the Sun whose quadrupole 
moment is only J20 < 2.3 x 10"'' [7, 8]. Moreover, in the case of a photon propagating near the planet the interaction 
between the gravitational field and the photon can no longer be considered static, because the planet moves around 
the Sun as the photon traverses through the Solar system [9, 10]. The optical interferometer designed for the space 
astrometric mission SIM [11] is capable of observing optical sources fairly close in the sky projection to planetary 
limbs with a microarcsecond accuracy. Similar resolution can be achieved for radio sources (quasars) with the Square 
Kilometer Array (SKA) [12] if it is included to the inter-continental baseline network of VLBI stations [13]. The 
Gaia [14] and OBSS [15] astrometric projects represent another alternative path to microarcsecond astrometry [16]. 
It is a challenge for the SIM and SKA interferometers as well as for Gaia and OBSS to measure the gravitational 
bending of light caused by various planetary multipoles and the orbital motion of the planets. This measurement, if 
sussessful, will be a cornerstone step in further deployment of theoretical principles of general relativity to fundamental 
astrometry and navigation at a new, exciting technological level. 

The first detection of gravitational bending of light by Jupiter was conducted in 1988 March 21 [17] (see also [18]), 
and the deflection term associated with the monopole field of Jupiter was determined to an accuracy of ~ 15% to be 
in agreement with Einstein's general relativity theory. Later on, the Hubble Space Telescope was used to measure the 
gravitational deflection of light of the bright star HD 148898 as it passed within a few seconds of arc near Jupiter's limb 
on 24 September 1995 [19]. Kopeikin [20] proposed to use Jupiter's orbital motion in order to measure the retardation 
effect in the time of transmission of the gravitational perturbation from the Jupiter to the photon, that appears as a 
small excess to the Shapiro time delay caused by the change in the gravitational force due to the motion of the planet. 
This proposal was executed experimentally in 2002 September 8, and the retardation of gravity (as compared with the 
time of propagation of light from the Jupiter to the observer) was measured to ~ 20% accuracy [21]. The retardation 
of gravity is determined by the ultimate value of the speed of gravity that was found to be equal to the speed of light 
as implied by general relativity. Physical interpretation of the gravity retardation effect crucially depends on the way 
of understanding of the theoretical concepts of general relativity and on how these concepts and principles are applied 
in experimental gravitational physics. Different authors elaborate on the multi-faceted realizations of these principles 
and depending on the approach tend to emphasize various physical aspects of the jovian experiment [22-29]. 

The jovian experiment of 2002 stimulated researchers from the Gaia consortsium [14] to take a next step in exploring 
the gravitational bending of light by major planets. Most notably, Crosta & Mignard [30] proposed to measure 
the deflection-of-light term associated with the axisymmetric (quadrupolar) part of Jupiter's gravitational field [31, 
32]. Their work is aimed at converting the earlier theoretical calculations [33-35] of light bending by gravitational 
multipoles into a practical algorithm for Gaia, thus, extending the relativistic techniques of astrometric data reduction 
worked out in a number of previous papers [9, 10, 36-39]. Detection and precise measurement of the quadrupolar 
deflection of light in the Solar system is important for providing an independent experimental support for the theory 
of gravitational lensing by clusters of galaxies. The cluster's gravitational potential (including invisible dark matter) is 
reconstructed from the observed distortion of images of background quasars under the assumption that the multipolar 
field of the gravitational lens deflects light exactly as predcited by general relativity [40, 41]. We emphasize that this 
assumption is highly plausible but still pending experimental confirmation that might be crucial for getting an unbiased 
estimate of the amount of dark matter in the universe. 

The work [30] is stimulating but it is incomplete and should be extended in several respects. First, it assumes that 
light propagates in the field of a static planet while the Jupiter moves on its orbit as light traverses the Solar system 
toward the observer. The authors of [30] do not specify at what particular instant of time the planet is fixed on its orbit 
and whether the orbital velocity of the Jupiter is essential for proper interpretation of the microarcsecond observation 
of light deflection. Second, Crosta & Mignard [30] implicitly assume that the center of mass of the planet deflecting 
light rays coincides precisely with the origin of the inertial coordinate system in the sky used for interpretation 
of the apparent displacements from the gravity-unperturbed (catalogue) positions of stars. This makes the dipole 
moment, P, of the gravitational field of Jupiter vanish, which significantly simplifies the theoretical calculation of 
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light bending. However, the assumption of J' = is not realistic because the instantaneous position of the planet's 
center of mass on its orbit is known with some error due to the finite precision of the jovian ephemeris limited to a 
few hundred kilometers [8]. The ephemeris error will unavoidably bring about a non-zero dipole moment that must 
be included in the multipolar expansion of the gravitational field of the planet along with its mass and the quadrupole 
moment. In other words, any realistic set of measurements can only be adequately interpreted within a certain model 
of the relativistic deflection of light, which includes parameters accounting for a possible shift of the true position 
of the planetary center of mass from the assumed origin of the deflection pattern, and their relative motion. It is 
important to realize that introduction of a static local coordinate system linked to the global star reference frame 
at few microarcsecond level, will give rise to non-physical deflection patterns caused by the planetary motion, which 
must be sorted out in any realistic experimental setup. 

The dipolar anisotropy in the light-ray deflection pattern is a spurious, coordinate-dependent effect and, hence, 
should be properly evaluated and supressed as much as possible by fitting the origin of the coordinate system used for 
data analysis to the center of mass of the planet. Until the effect of the gravitational dipole is properly removed from 
observations it will forge a model-dependent quadrupolar deflection of light because of the translational change in the 
planetary moments of inertia - the effect known as the parallel-axis or Steiner's theorem [42] (see Eq. (17) below). 
This translationally-induced quadrupolar distortion of the light-ray deflection pattern should be clearly discerned 
from that caused by the physical quadrupole moment of the planet J2. This paper discusses the theoretical and 
observational aspects of monopolar, dipolar, and quadrupolar light-ray deflections. We investigate how the spurious 
deflections can be separated from the physical ones caused by the intrinsic quadrupole moment of the planet J2, and 
how to use the measurement of the dipolar anisotropy to test the relativistic effects caused by the time-dependent 
component of the gravitational field deflecting light rays. 

We describe our mathematical notations in section II. Propagation of light ray in flat space-time and parametriza- 
tion of the gravity-imperturbed light-ray trajectory is given in section III. Analytic description of the retarded 
gravitational field of a moving planet and its multipolar decomposition is introduced in section IV. Null geodesic 
equations describing the gravitationally-perturbed light-ray propagation, mathematical technique of their integration, 
and definition of the observed angle of the gravitational light-ray deflection are included to section V. Specific light-ray 
deflection patterns caused by mass, dipole, and quadrupole moments of gravitational field are outlined in section VI. 
We discuss observable relativistic effects in section VII, and technological capabilities of SIM and SKA interferometers 
to pursue astrometric measurements of these effects with microarcsecond accuracy, in section VIII. 

II. NOTATIONS 

To simplify our theoretical equations we use the system of units in which the fundamental speed c and the universal 
gravitational constant G are both equal to unity: G = c = 1. If necessary, these fundamental constants can be 
restored in equations by making use of dimensional arguments [43, 44]. 

We further assume that general relativity is valid, which implies that gravity operates on the null cone and the force 
of gravitational interaction propagates with the fundamental speed c [45] . We suggest that each photon incoming to 
the Solar system propagates in vacuum, and its physical speed is equal to c. It simplifies our consideration as we 
avoid discussing the perturbations in propagation of light due to dispersive effects in plasma. 

Latin indices take values 1,2,3, and the Greek ones run from to 3. The Kroneker symbol (a diagonal unit matrix) 
is denoted 6ij = diag(l, 1, 1), and the fully anti-symmetric symbol of Levi-Civita e^fc is defined in such a way that 
ei23 = +1- The Minkowski metric is rja^ = diag(— 1, 1, 1, 1). Greek indices are raised and lowered with the Minkowski 
metric rjajs- By convention, Latin indices are raised and lowered with the Kroneker symbol Sij which makes no 
difference between super- and sub-script Latin indices. Repeated indices indicate the Einstein summation rule, for 
example, rja^A^Bf^ = A^Bl^ = AoB° + A^B^ + AjS^ + A^B^. 

An arrow above a letter denotes a spatial vector, for instance x = x' = {x^ ,x^). A dot or a cross between 
two spatial vectors denote the Euclidean scalar or vector products respectively: A ■ B = AiB^ = AiBi = A^B\ and 
A X B = {A X By = CijkA^B^. Angular brackets around a pair of the Latin indices of a spatial tensor of rank two 
denote its symmetric and trace-free (STF) part, for example, [46] 

j<ij> = i ^jij ^ _ Igijjkk ^ 
2 3 

where = 5'=^/'=^. 

Partial derivatives with respect to four-dimentional coordinates x" are denoted with comma so that for any differen- 
tiable function F{t,x) one has F^^ = dF/dx". Partial derivatives of F{t,x) with respect to spatial coordinates x^ are 
denoted as Fj = dF/dx^ = VF. Partial derivative of function F{t,x) with respect to time is denoted Fq = dF/dt. 
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Total derivative of F{t, x) with respect to time is denoted with overdot, that is 



• dF dF dx 

F = — = 1 VF . (2) 

dt dt dt ^ ' 

Notice that, in general, F ^ dF/dt except when F is a function of time only. 

III. THE UNPERTURBED LIGHT-RAY TRAJECTORY 

Let us introduce a global coordinate system x" = (x*^, x*) = {t, x*) which coincides with an inertial frame at infinity 
where the gravitational field is absent. We shall assume that the coordinate system x" is static with respect to the 
observer but we shall prove in section V that our results are gauge-independent and Lorentz-invariant and, in fact, 
are valid in any other frame moving with respect to the observer with constant velocity. 

Let us consider a bundle of light-rays emitted by a source of light (star, quasar) simultaneously and moving as a 
narrow beam along parallel lines toward the Solar system. In the absence of gravitational field each light particle 
(photon) from the bundle propagates in the coordinate system x" along a straight line 

x' =xi + k\t-to) , (3) 

where to and Xq = x^{to) are the time and space coordinates of the photon at the time of emission, and fc' is the 

unit vector along the unperturbed photon's trajectory (see Fig. 1). We assume that the photon hits the detector (is 
observed) at time ti when its coordinate x\ = x^{t\). Let us denote the time of the closest approach of the photon to 
the origin of the coordinate system as [47] 

= ^0 — • xo . (4) 
It is mathematically convenient [34] to introduce a parameter r along the unperturbed light ray 

T = t — tt, , (5) 

and a constant (light-ray impact-parameter) vector 

f = piix^ = piixl ' (6) 

that points out from the origin of the coordinate system to the point of the closest approach of the unperturbed light 
ray (see Fig. 1). Here P'^ = (5'^ — fc'fc-' is the operator of projection onto the plane of the sky, that is orthogonal to 
the vector fc*. By definition, P^i pi^ = P'*^. Parametrization of the unperturbed light-ray trajectory given by Eqs. 
(4)-(6) converts Eq. (3) to 

= kw + e , (7) 

where t and are independent of each other, and can be considered as coordinates in the 2+1 manifold of the 
light-ray bundle because the projection (6) makes have only two independent components. 

At each instant of time t, the distance r of the photon from the origin of the coordinate system is 



r = + rf2 , (8) 

where d = \^\ is the absolute value of the impact parameter of the photon, that is constant for each light ray from the 
parallel beam of the light-ray bundle moving in the direction k. Notice that definition (5) implies that r = when 
t = t*, r < when t <tt:, and positive otherwise (see Fig. 1). 



IV. THE GRAVITATIONAL FIELD 



A. The field equations 

We shall consider the Solar system to be isolated and space-time to be asymptotically flat which means there are 

no other masses outside of the Solar system. The gravitational field of the Solar system is produced by the Sun and 
the planets which curve space-time and deflect light. In what follows, for the sake of simplicity we are taking into 
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account the gravitational field of two bodies only - the Sun and a planet moving around the barycenter of the Solar 
system. Futhcrmore, wc consider the solar gravitational field as spherically-symmetric in its own, proper reference 
frame because in this frame the dipole and quadrupole moments of the Sun are negligible. Gravitational deflection 
of light by the Sun is well-known [48] and we will basically focus on the discussion of the light bending by a moving 
axisymmetric planet. 

The planet moves around the barycenter of the system as a light ray propagates from a star toward the observer. 
Position of the planet with respect to the origin of the coordinates at time t is defined by vector xp = xp{t), 
its velocity is denoted vp = dxp/dt, and acceleration ap = dvp/dt. We assume that the planet is axisymmetric 
around the unit vector s = s{t) that defines the rotational axis of the planet at time t. This vector can change its 
orientation in space due to precession. Our calculation method is general enough, and we do not need to assume that 
the parameters characterizing translational and rotational motion of the light-ray deflecting bodies are either constant 
or equal to zero, as it was postulated in [30] where the authors assumed that xp = 0, along with vp = ap = and s 
is a constant vector. We shall calculate relativistic deflection of light by a planetary axisymmetric gravitational field 
but ignore the relativistic effects caused by the gravimagnetic field due to the intrinsic rotation of planet, since it is 
negligibly small and can not be detected at the microarcsecond resolution [10]. We refer the reader to paper [5] where 
detailed estimates of the magnitude of the gravitational light deflection caused by various parameters of the Solar 
system bodies are given in Table 1. 

The gravitational field of the Solar system is described by the metric tensor 

9al3 = Vaff + ha0 , (9) 

where rjafj is the constant Minkowski metric, and hafj = hafj(t,x) is its perturbation which is associated in general 
relativity with gravitational potentials. Let us impose the harmonic gauge condition on the potentials 

- \ri"^Kli = , (10) 

which is convenient and mathematically powerful choice for solving Einstein's gravity field equations [49, 50]. Outside 
of the planet, in vacuum, and in the harmonic gauge (10), the linearized Einstein equations for the field ha0 are 
homogeneous wave equation [46, 51] 

{-^+^')h„0 = O. (11) 

A question can arise about universality of the harmonic gauge and physical interpretation of the wave solution of the 
gravity-field equation (11). Indeed, Einstein's theory of general relativity is formulated in a covariant tensor form while 
the wave equation (11) is valid in a particular harmonic gauge only. First of all, we notice that the harmonic gauge 
is not reduced to a single coordinate system but admits a whole class of both global and local harmonic coordinates 
related to each other by coordinate transformations which do not violate the harmonic gauge condition (10). The 
class of the global coordinates consists of the reference frames which are moving with respect to each other with 
constant velocities in the asymptotically flat space-time. These asymptotically-inertial reference frames are connected 
thorugh the Lorentz transformation as was shown by Fock [50] for a generic case of gravitational field including 
all non-linearities of Einstein's equations. Local harmonic coordinates were introduced to gravitational physics by 
Thorne and Hartle [52] . The law of transformation between the local harmonic coordinates extends the linear Lorentz 
transformation to the class of harmonic polynomials that arc solutions of the homogeneous wave equation [53 58]. 
They have a great practical value for modern fundamental astronomy [59, 60] . Second important observation is that 
the class of the harmonic coordinates is actually much more generic than physicists used to think. In particular, we 
have shown [2, 61, 62] that the standard Arnowitt-Deser-Misner gauge [63] widely used for dynamic formulation of 
general relativity, may be viewed as a sub-class of the harmonic gauge suitable for description of gravitational field 
of an isolated astronomical system in vacuum. Final remark is that we always work in this paper with solutions 
of the gravity-field wave equation (11) in terms of an observable quantity that is the deflection angle of a light ray 
referred to asymptotically-flat space-time. This quantity is invariant with respect to the gauge transformations and, 
hence, all our results obtained in this paper are valid in arbitrary gauge. The harmonic gauge is simply a convenient 
mathematical tool facilitating calculations and physical interpretation of the observed relativistic effects. 

B. The gravitational multipoles 

A general solution of Eq. (11) is given in the form of a multipolar expansion [46, 51] depending on time-dependent 
parameters characterizing the intrinsic multipole moments of the planet as well as on the displacement xp of its center 
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of mass from the origin of the coordinate system. We shall omit in this paper the relativistic effects of the planet's spin, 
and consider only the mass monopole, dipolc, and quadrupole terms in the miiltipolar expansion of the gravitational 
field. Relativistic effects caused by the spin and higher-order mass multipoles are discussed in [10, 34, 61, 62]. 

Because Eq. (11) is linear, we can consider the gravitational field of the Solar system to be a linear superposition 
of the individual fields. For each massive body the solution of Eq. (11) is [46, 51] 



2M 
r 

4P{s) 



d 



2r{s) 



+ 



+ 



/<*J>(s) 



(12) 

(13) 
(14) 



where r is the distance from the origin of the coordinate system to the field point x" = (t,.T*), and the angular 
brackets around two indices denote the STF tensor as explained in Eq. (1). Notice that the second time derivative of 
the dipole moment P is absent due to the law of conservation of the linear momentum of the body [46, 51]. 

The gravitational field described by Eqs. (12)-(14) should be interpreted as that created by a massive planet placed 
at the origin of the coordinate system and characterized by several parameters, of which M is a constant mass of the 
planet, P{s) is the dipole, and /^*^^(s) is the quadrupole moment - both taken at the retarded instant of time 



s = t — r . 



(15) 



The retardation is a direct consequence of the retarded (causal) solution of the gravitational wave equation (11). The 
null cone corresponding to the causal region of influence of the gravitational field of a planet is shown in Fig. 2. The 
planet moves along its own world line so that observer is expected to see any effect of the planet's gravitational field 
with the retardation that is due to the finite speed of the propagation of gravity on the future part of the null cone. 
Figure 2 gives a general-relativistic picture of the process of propagation of gravity from the moving planet to the 
observer but it is also valid in a number of alternative theories of gravity [22, 64]. 

The dipole and quadrupole moments are fully taken into account because they produce a relativistic defiection of 
light that can be of significant value for several planets [5, 30]. In general relativity, all relativistic efi'ects due to 
the dipole moment of the gravitational field are coordinate-dependent and, hence, can be eliminated if the origin 
of the coordinate system is placed exactly at the center of mass of the real planet. However, we retain the dipole 
moment /* in the multipolar decomposition of the planetary gravitational field (12)-(14) because in practice we do not 
know accurately enough where the center of mass or the planet is and, hence, one needs to parametrize its position. 
Furthermore, we should consider that even if we were able to to render = at a particular instant of time, the 



planet has a residual velocity with respect to the coordinate frame a;", and its center of mass. 



is not fixed at 



the origin of this frame for any observation of finite duration. It is therefore, important to properly parametrize and 
evaluate the infiuence of the dipole moment on the light deflection. 

If the center of mass of the planet is shifted from the origin of the coordinate system by a spatial distance L' = Xp, 
the dipole and quadrupole moments of the gravitational field entering the metric tensor perturbations (12)-(14) are 
defined in the linearized approximation of general relativity by the following equations 



P = 

J<i3> _ 



Mx'p , 

J<^3> + 



Mxi'x'ji' 



(16) 
(17) 



where J^'^^ is the intrinsic quadrupole moment of the planet in its own proper reference frame whose origin coincides 
with the planetary center of mass, and we made use of the parallel-axis theorem [42] to split to J'^ and the 



translational term 



Mxf 



Equations (16), (17) approximate more precise, post-Newtonian definitions of the 



multipole moments used in the relativistic celestial mechanics of N-body system [58] . However, the post-Newtonian 
corrections to equations (16), (17) yield terms in the light deflection which are much less than 1 /uas, and hence, can 
be ignored in the present paper. 



Both the dipole, Mxp, and the translational quadrupole, 



Mx<' 



will affect the post-fit results of our light-ray 
deflection measurement, and thus, are dircc;tly observable in the light-ray deflection pattern. But these deflections are 

spurious, coordinate-dependent effects and the goal of the data analysis system is to reduce their impact on the values 
of the fitting parameters within the accuracy of astrometric observations. To this end the dipolar and quadrupolar 
deflections associated with the translation = Xp must be incorporated to the data analysis in order to assume full 
control on their influence on the values of the fitting physical parameters like J2, etc. If this is not done, the coordinate- 
dependent effects can not be eliminated and the light deflection due to the dipole, Mxp, and the translational 
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quadrupolc, Mx^'^x-'p , may exceed the physical deflection of light caused by the intrinsic quadrupole moment, J2, of 
the planet, thus, making its measurement in the gravitational light-ray deflection experiments impossible. 

In what follows, we assume that the planet is axisymmetric around its rotational axis defined at each instant of 
time by a unit vector s*. The planet has equal equatorial moments of inertia A = B, and the axial moment of inertia 
C ^ A. The dynamic oblateness of the planet is denoted as J2 = (C — A) /A [48, 65]. This definition yields the 
intrinsic quadrupole moment represented as an STF tensor of the second rank [65] 

J<^j> = MJ2R^ (^s's^ - y^^^ , (18) 

where R is the equatorial radius of the planet. It is immediately seen from Eq. (18) that the quadrupole moment 
is symmetric and trace-free, that is J<"^ = 0, in accordance with its definition. We also notice that the first time 
derivative of the dipole and quadrupole moments 

P = Mv'p , (19) 
7<'^> = 2Mx<'vj^ , (20) 

should be used in calculations of hoi in Eq. (13). We shall neglect the first time derivative of J^'-'^ in hoi, and 
the second time derivative of the overall quadrupole moment in hij as they lead to higher-order relativistic 

effects in light deflection. For example, as follows from the subsequent calculations, the light-ray deflections caused by 
the second derivative of the quadrupole moment have magnitude of the order of au^vplc)^ , where au is the 

gravitational bending due to the mass monopole of the planet. In case of Jupiter the maximal value of olm ~ 16300 
/uas [5, 30], and vp/c ~ 4.5 x 10~^ [20], hence, cxMivpjc)'^ <C 1 /(xas. 



V. THE LIGHT-RAY GRAVITATIONAL PERTURBATIONS 



A. The light-ray propagation equation 
The equation of motion of light in gravitational field is given by the light-ray geodesic [43, 44] 

-^^Vl^K^^K^ = (21) 
where = dx°' /d\ is the wave vector, A is the affine parameter along the light ray, and 

= la"^ {90^,u + 9P^,^ - 5^.,/3) . (22) 

is the Christoffel symbol [43, 44]. 

When one substitutes the expansion of the metric tensor (9) to Eq. (21), and transforms the affine parameter A in 
this equation to the coordinate time t, the light geodesic equation (21) is reduced to the following three-dimensional 
form [48] 



^2 i 

F\t,x), (23) 



df^ 

where the gravitational force [48] 

-F* = -hoo,i — hoifi — -^hoo^ox^ — hik^x'^ — (ft-oi^fe — hok,i)x^ — (24) 

hoQ^hX^x' - (^hikj - ^hkj,ij x'^x^ + {^hkjfi - hok,j^ x'^x^x' , 

depends on the perturbation of the metric tensor hajs and the coordinate velocity of photon x*. We notice that making 
use of time t in place of the parameter A does not change direction of propagation of the light ray, and is merely 
a technical tool that allows us to express coordinate of the photon as a function of the same time argument which 
governs evolution of the multipole moments of the gravitational field. 

Equation (23) has an unperturbed {F^ = 0) solution described in section III. In a weak-gravitational field approx- 
imation the light-ray geodesic equation (23) has a unique solution given by 

a;«=fcV + f + S'(t,0, (25) 
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where is a small perturbation of the unperturbed light ray trajectory. In terms of the parameter r and the 
projection coordinates equation (23) is significantly simplified, so that the perturbation obeys the following 
ordinary differential equation [2, 61] 

^ = \k'^kP^-^ - ^ [k-h,^ + \k^h,, - ^k^k^kPh,,^ , (26) 

where k" = (l,fc*) is a null vector along the unperturbed trajectory of the light ray. Equation (26) describes two 
relativistic effects - gravitational bending and time delay of light which are closely related to each other [62] . In this 
paper we discuss only the deflection of light. 



B. The null cone integration technique 



Without any restriction on the ratio of the impact parameter d of the light-ray trajectory to the distance ri from 
the origin of the cioordinatc system to observer, the total angle of the gravitational deflection of light in the plane of 
the sky is given by a vector [2, 36] 



a' 



d"i 
dr 



(27) 



where Aa' includes relativistic corrections due to both the finite distance ro of the source of light from the planet, the 

finite distance ri between the observer and the planet, and the difference of the local inertial frame of observer [66, 67] 
from the coordinate system = {t,x^) introduced for the calculation of light-ray propagation. We further assume 
that ro ^ oo, and the observer is at rest at a sufficiently large distance from the light-ray deflecting body which are 
realistic assumptions [68, 69]. Under these circumstances Aa* is negligible, and can be omited. Indeed, the largest 
contribution to Aa* associated with the flnite distance ro of the star from the Solar system (the planet), is smaller 
than the first term in the left side of equation (27) by a factor of ri/ro [2, 9] that, in case of the Jupiter, is about 
2 X 10~^ even for the closest star a Centauri, making Aa' <C 1 /txas [70]. Contribution to Aa' due to the difference 
between the local inertial and global coordinate frames is smaller than the first term in the left side of equation (27) 
by a factor of {d/ri){vp/cf' [2, 9] that also amounts to jAa^j ^ 1 /ias, and can be ignored. In cases when Aa* may be 
important, for example, in studying refraction of light by gravitational waves [2, 62] or in cosmological gravitational 
lensing where ro 7^ and ri is comparable with ro, this term is solved in our papers [9, 10, 61]. 

Integrating Eq. (26) along the unperturbed light ray trajectory and substituting the result to Eq. (27) yields [9, 10] 



1 d f^^ -> r 

a* = "2^ J k"k^h„/}iT,^)dT + k"'P'^^hj„{si,xi) - hj„iso,xo) 

where the integrand is taken on the unpertubed light-ray trajectory, and can be represented as 



2P{s) 



+ 



r 

I<ip>{s) 



d_ 
' dT 

d^ 
dr^ 



2kU^{s) 



(28) 



(29) 



2ykPI<jP>{s) 



d d 

dT 



2fcP7<jp>(s) 



the limits of integration ti and tq are the values of the parameter t taken at the time of observation, ii, and emission, 
fo, respectively (notice that To < and ri > 0), and s = + r — r with r = r{T, ^ ) given by Eq. (8). Integration in the 
right-hand side of Eq. (28) can be performed easily if one adopts that in general relativity the speed of propagation 
of gravity and the speed of light are numerically the same, making the light cone hypersurface be coinciding with one 
of the gravity null cone, so that in practical experiments gravity interacts with light in a way shown and explained 
in Fig. 3. In some of the alternative theories of gravity the hypersurfaces of the light and gravity null cones do not 
coincide as the speed of propagation of gravity and light are different [25, 71-74]. We have discussed the gravitational 
deflection of light caused by moving bodies, in such theories in a number of other papers [22, 64, 72] which can be 
used as references. We do not discuss the deflection of light in the alternative theories of gravity in the present paper. 

We have found [2, 61, 62] that for any smooth function of the retarded time, F(s)/r, where r = ^t'^ + d"^ and 
s = t — r = U+ T— Vt^ + d'^, the integration can be performed explicitly, and it yields 



Tn 



_d_ 

' A. 



Fis) 



Fis) 



dT 



dT = 



d_ 
n 



Fis) 



dT- 



(1 cosx)-F(si) - (1 - cosxo)F(so) 



d? 



ro 



(30) 
(31) 
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where the retarded times 



So = to 
si = h 



(32) 
(33) 



cosxo = I'ToI/t'o, cosx = |Ti|/ri arc referred to distances ro = |x(io)|, fi = |^(ti)| of the photon from the origin of the 
coordinate system taken at the times of emission and observation of Ught respectively (see Fig. 1). Signs in the right 
side of Ecj. (30) are valid in case shown in Fig. 1 when the light-ray deflecting body is located between the source of 
light and observer. If the observer is located between the source of light and the light-ray deflecting body, Eq. (30) 
is still applicable after replacing cos x ^ " cos x [2] ■ Finally, if d = Vi, that is the light ray arrives from the direction 
perpendicular to that to the light-ray deflecting body, we must take cosx = in Eq. (30). All these situations are 
carefully discussed in [2] to which the reader is referred for a comprehensive theoretical review. In this paper we 
shall discuss the configuration shown in Fig. 1 because only this case is practically important for the gravitational 
deflection of starlight by the Solar system planets. 

It is noted that the retarded time Eq. (33) describes the null direction connecting the planet and the observer and 
lying on the future gravity null cone with the planet at its vertex as shown in Fig. 3. Light moves along a different 
null direction connecting the star and the observer and lying on the past light cone of the observer. Therefore. Eq. 
(33) describes the retardation effect in propagation of gravity force from the moving planet to a photon as the photon 
propagates toward the observer and is subsequently detected at time ti . This retarded component in the interaction 
of gravity with light was measured within 20% in the jovian VLBI experiment on September 8, 2002, and led to the 
direct observational confirmation of the general relativistic postulate that the speed of gravity and light are the same 
within the observational error [20-22]. 



Review article [25] argues that the speed of propagation of the gravity force from the planet to the traversing 
photon is irrelevant in the light-ray deflection experiments done in the Solar system with a uniformly moving planet, 
and that the light deflection angle can be calculated correctly even if the speed of gravity is infinite. In other words, 
the author of [25] mantains the point that the causal structure of the future null cone of the gravity field (see Fig. 
3) is not essential to the calculation of the light-ray deflection angle so far as acceleration of the planet is ignored. 
This point of view stems from the belief, explicitly formulated and adopted in the PPN formalism [73], according to 
which the principle of equivalence demands nothing about the speed of gravity that determines the causal space-time 
hypersurfacc on which gravity propagates. We are told that this is because the principle of equivalence operates with 
derivatives of the metric tensor of the flrst order (the Christoffel symbols) while the speed of gravity is determined 
only by the second time derivatives of the metric that appear in the fleld equations of each metric theory of gravity 
(for more detail see section 10.1 in [73]). 

This PPN postulate disentangles gravitational field from the geometry of the space-time manifold and ignores 
existence of other (global or geometric) techniques [43] to study the causality of gravitational field which is a direct 
consequence of its ultimate speed of propagation. The causal structure of gravitational field is derived not so much 
from the field equations but from exploration of behavior of a set of future directed timelike and/or null geodesies 
in a given space-time manifold [75]. This behavior can be determined, at least in a close neigborhood of any event, 
by the geodesic equations without imposing Einstein's field equations [76] . The null geodesies define the causal past 
of observer, that is the region bounded by the past light cone in Fig. 3, and the causal future of the gravitational 
field of a massive body, that is the region inside the future gravity cone in Fig. 3. General- relativistic gravitational 
interaction of photon with the body implies that the causal past of the observer must touch the causal future of the 
gravitational field along a null direction as demonstrated in Fig. 3 [77] . 

The equivalence principle tells us that in a local reference frame light moves along a straight line [43, 44]. This 
implies that in the global reference frame the light-ray trajectory is bent because the local reference frame falls in the 
gravitational field of a global frame with acceleration, and the space is curved [43]. Testing gravitational light-ray 
deflection is not an experiment in a single local frame of reference as photon propagates through a continuous sequence 
of such local frames, thus, accumulating the pointwise influence of the gravitational fleld at different parts of the light- 
ray trajectory. Photon's wave-vector propagator is an integral of the affine connection that is the gravity force 
shown in the right-hand side of Eq. (23). This propagator contracted with the projection operator P^^ = 5^^ — k^k^ , 
yields the integrated deflection angle a*. The magnitude and/or direction of the deflection angle are functions not only 
of the wave vector of observed starlight, but of the position of the light-ray deflecting body with respect to observer 
as well. If one is able to derive position of the massive body from the precise measurement of the angle a* and to 
confirm that the body and the observer are connected by a null line which is a characteristic of the future gravity null 



C. The speed of gravity, causality, and the principle of equivalence 
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cone, it gives a direct proof of the causal nature of the gravitational field and allows us to measure its finite speed 
of propagation [20]. However, this measurement of the causal property of gravity is impossible if the gravitational 
field is static. Indeed, in case of a static planet the gravitational field does not change as photon moves from star to 
observer, and the causal character of the gravitational field gets hidden from observer because the planet is always at 
a fixed distance from the observer (see Fig. 4a). 

The situation changes dramatically if the planet moves with respect to observer because it makes the distance 
between the planet and the observer depending on time (see Fig. 4b). In this case, photon traverses through the 
gravitational field that changes on the light-ray trajectory due to the planetary motion, even if this motion is uniform. 
Had the speed of the propagation of the planet's gravitational field been different from the speed of light it would have 
unavoidably led to violation of the causal nature of gravity that would be inconsistent with the principle of equivalence 
[74] . Indeed, in this case the light null cone would be different from the gravity null cone, and the moving planet would 
not deflect the observed photon from the retarded position as predicted by general theory of relativity which predicts 
that the angle of the gravitational deflection of light is a gauge-invariant quantity (compare Fig. 4a with Fig. 4b 
which are connected via Lorentz transformation). For example, the instantaneous propagation of gravity would imply 
that one could determine current (as opposed to retarded) position of the planet on its orbit from observation of the 
gravitational deflection of light, that is the gravitational field would transmit information about the planet's spatial 
location to observer faster than light. This violates the principle of causality, and we conclude that correct description 
of the gravitational physics of the light-ray deflection experiment with a moving planet requires taking into account 
both the light and gravity null cones as demonstrated in Figs. 2, 3, 4, and supported by calculations in this paper as 
well as the discussion given in [72] . Additional arguments backing up the concept of the retardation of gravity in the 
light-ray deflection experiments by moving planets are discussed in our paper [22] both in the framework of general 
relativity and in a bi- metric theory of gravity proposed by Carlip [80]. The present paper and calculations given in 
[22] refute the superficial understanding of gravitational physics of the light-ray deflection experiments advocated in 
the PPN formalism [25] and in [80] . The principle of equivalence does imply the principle of causality for gravitational 
fleld and demand the ultimate speed of gravity to be equal to the speed of light which can be (and was) confirmed in 
the Solar system experiments of starlight deflection by a moving planet [21]. 



VI. LIGHT-RAY DEFLECTION PATTERNS 



A. The deflection angle 



In calculation of the light-ray deflection angle all terms related to the time of emission are negligibly small because 
1 — cosxo — '^^/(2?'o) ^ 1) ^iid can be omitted since we assume that ro OC'. We also assume that only the stars 
that are sufficiently close to the planetary limb in the plane of the sky are observed. This makes the angle x <C 1, 
and hence cosx — 1 — rf^/2ri and sinx = d/n <^ 1. We shall neglect all terms proportional to sinx and approximate 
cos X = 1 everywhere but in the calculation of the monopolar deflection, where more exact approximation of cos x is 
required. 

Integrating Eq. (28) with the technique shown in Eqs. (30)-(31), and keeping only the leading terms, we deduce 

a'=aQ + ail + + , (34) 

where a\f, a^, Uq are the angles of the relativistic deflection of light caused by the planetary mass monopole, dipole, 
and quadrupole moments respectively, and a© stands for the well-known deflection angle caused by the Sun's gravity 
[48] . The deflection terms are defined by the following equations 



a 



M 



2M(l + cosx)|2 ' 
-2fc^P(5i)(l + cosx)|^-2A 



(H-cosx) 



a, 



Q 



(1 -hcosx) 



^P/<^P>(si) 
d2 



(35) 
(36) 

(37) 



of which Eq. (35) is exact, and in Eqs. (36), (37) the terms proportional to sinx and higher, have been omitted. 
These residual terms are of the order of aod/ri and aQd/r\, that is negligible, since max[a£),Q;Q] ~ 2000 /uas (see 
Figs. 6, 7), and d/ri ~ 10~^ for Jupiter and about 10~^ for Saturn. 
Taking all the partial derivatives with accounting for 

d cos X C 



cosxsm 



d 1 
^,«-d 



(38) 
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and approximating cosx = 1, sinx = in the dipolar and the quadrupolar deflections, we reduce Eqs. (35)-(37) to 

AM 

= — (l+cosx)n% (39) 

otjy = — -^-^ [vH-n^ - wH'm?) — -n^ , (40) 

Uq = — - (n^'n^n^ - n^'m^mF - m''m^n'' - m^mFn') , (41) 



where P and are determined by Eqs. (16) and (17) respectively, the unit vectors n* = ^V*^; = {k x fiy, and 
again all terms proportional to sin x have been neglected. On the other hand, the term with cos x in Eq. (39) must 
be retained as it reach magnitude of the order of 2 /xas when x = 90° for Jupiter, and is about 1 /uas when x = 16° 
for Saturn, that is the small-angle approximation in Eq. (39) is invalid. 

It is evident that the deflection angles a)j and ag depend on the value of the dipole and quadrupole moments of 
the planet taken at the retarded time si ~ ti — ri, where ti is the time of observation. The deflections caused by the 
translation x^p are unphysical and can be removed, at least in theory, by chosing the origin of the coordinate system 
used for calculation of the deflection angles at the center of mass of the planet taken at the retarded instant of time 
with respect to the observer. This situation is shown in Figs. 3, 4b where the past light cone is constructed from the 
null characteristics of the light rays coming from all the stars to the observer, and the future gravity null cone is made 
of the null characteristics of the gravity-field wave equation (11) that demarcate the region of the causal interaction 
of the gravity force of the moving planet with any other particles inchiding photons. 

The deflection angle a' is an observable and gauge-invariant quantity measured with respect to the unperturbed 
direction fc' of incoming photon defined in the asyniptotically-fiat space-time (at past null infinity). For this reason, 
the relativistic efi'ects associated with the gravitational deflection of light have direct physical interpretation that is 
discussed in the rest of the present paper. 

B. Snapshot patterns 

Substituting Eqs. (16)-(19) to Eqs. (40), (41) yields 

O.M = a(l+cosx)n, (42) 
L 



an = a — 
a 



{z- n)n — [z- rh)rh , (43) 

i?2 



dq = aJ2-^| [(s • n)^ - (s • to)^] n - 2(s • n)(s • TO)m| (44) 

rf2 



-|- a [{z- nf — {z- m)2] n — 2{z- n){z- m)mj , 



where 



d 

is the light-ray deflection angle caused by the planetary mass, and 



a = aiimb-T , (45) 



*^limb 



A{l-k-vp)^ (46) 



is the deflection angle on the planetary limb, L = \xp{si)\ is the absolute value of the displacement of the planetary 
center of mass from the origin of the coordinate system, the unit vector z = xp{si)/L points from the origin of the 
coordinate system toward the center of mass of the planet, and all time-dependent quantities like xp, z, and s, are 
computed at the retarded time Si = ti — ri corresponding to the time of observation ti. 

It is possible to prove [81] that Eqs. (42)-(44) are still applicable in the class of scalar-tensor theories of gravity 
after a formal replacement: a ^ = (l + 7)01/2, where 7 is the PPN parameter characterizing the impact of the 
scalar field on the light propagation [25]. We have also moved the term depending on the first time derivative of 
the dipole moment from Eq. (43) to Eq. (46). This velocity-dependent term is purely radial and corresponds to 
the relativistic correction of the planetary mass M caused by the radial Doppler shift. This gravimagnetic correction 
correlates with the PPN parameter 7 setting a natural limit on the precision of its measurement from a single-epoch 
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experiment. In general relativity, 7=1, and the maximal general-relativistic deflection of light on the planetary limbs 
are respectively: aumb — 16280 /ias for Jupiter, and aiimb — 5772 iias for Saturn [5, 30]. Taking into account in Eq. 
(46) the gravimagnetic correction to the total mass makes these deflections modulated up to a few /uas, depending on 
the relative orbital velocity of the planet and the observer. 

It is worthwhile to emphasize that the entire theory of the microarcsecond light-ray deflection given in this paper is 
gauge- independent and Lorentz- invariant, that is all equations are valid in any arbitrary frame moving with respect 
to the original coordinate system = {t,x^) with velocity /3\ Equations (42)-(44) of the light-ray deflection angles 
remain invariant with the velocity replaced to — Lorentz-invariance of Eqs. (42)-(44) can be easily understood 
if one remembers that all quantities entering these equations are referred to the events connected by the null cones 
shown in Figs. 2, 3 that are, of course, invariant with respect to the Lorentz transformations. 

Eqs. (42)-(44) define the deflection patterns caused respectively by the monopolar, dipolar, and quadrupolar 
components of the planetary gravitational field with respect to the coordinate system a;" = {t,x'^). Because velocity 
of motion of planets in the Solar system are much smaller than the speed of light/gravity we can assume that during 
the time of propagation of photon across the Solar system the maximal displacement Xp of the center of mass of the 
planet that is allowed by the ephemeris error or some other violation in the data processing computational approach, 
is smaller than the impact parameter, L/d < 1. Let us consider the light-ray deflection patterns for a single epoch of 
observation, ti , mapping the deflections in the plane of the sky around the planet like it will appear in a photographic 
snapshot. To vizualize the instantaneous two-dimensional deflection patterns, it is instructive to project Eqs. (42)- 
(44) on the directions that are parallel and orthogonal to the ecliptic plane of the Solar system as shown in Fig. 5. 
The unit vectors n' and m' are expanded in such local tangential coordinates into 

n = Bfj cos ip + e\sunp , (47) 
m = —6/3 sin -|- ex cosy , (48) 

where 1^ is called the position angle [65], the unit vectors e\ and e}? are mutually orthogonal and directed along the 
increasing ecliptic longitude A and latitude (3 respectively. The unit vectors z and s are expanded in this tangential 
frame as follows 

z = ep cos + eA sin ^ , (49) 
3 = 6/3 cos u> + exsinijj , (50) 

where ip and oj are the position angles of the unit vectors z and s respectively. 
Substituting Eqs. (47)-(50) to Eqs. (42)-(44) yields 

cJm = q(1 + cosx) (e/3 cos</? + ca sin(y5) , (51) 



L 

ao = a — 

a 



Cfj cos{'ip — 2ip) — e\ sin(^ — 2ip) , (52) 



i?2 



dq = OiJ"!-^ cos(2'^ - 3(^) - e\ sin(2V' - 3(p) (53) 

6/3 cos(2cj — 3</?) — e\ sin(2u' — 3</?) 



In this form, the equations clearly show that the directions of deflection angles du, d]j and dg vary proportional to 
simp, sin2ip and sin3(^, respectively, as one goes around the planetary center. The amplitudes of q;m, c«d, and ag 
fall off as 1/x, 1/x^, and l/x^, respectively, where x is the impact angle in the plane of the sky between the star and 
the coordinate system origin (see Fig. 1). 

Snapshot patterns of deflections du, do, and dq for a stellar field surrounding Jupiter and Saturn are shown in 
Figs. 6, 7 where we use color grades to denote the magnitude of the deflection as a fmiction of the impact angle x of the 
light ray, and arrows to specify the direction of deflection, which depends on the position angle in accordance with 
Eqs. (51)-(53). The patterns are shown for three difi'erent magnitudes L of the displacement vector xp: L\ = 350 km 
(left), L2 = 3500 km (middle), L3 = 35000 km (right) for Jupiter, and Li = 250 km (left), L2 = 2500 km (middle), 
and L3 = 25000 km (right) for Saturn. 

The monopolar deflection of light cSm is depicted in subplots (a) of Figs. 6, 7. It is purely radial (that is aligned with 
vector n), and does not depend on the magnitude and direction of the displacement vector xp. The dipolar deflection 
of light, do, is directly proportional to the magnitude L of the displacement, and its orientation is determined 
by vector xp in the plane of the sky. The quadrupolar deflection of light, dq, is quadratically proportional to the 
magnitude L of the displacement between the origin of the coordinate system and the planet's center of mass as well as 
on the magnitude of the intrinsic quadrupole moment of the planet J2. The orientation of the quadrupolar deflection 
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pattern is determined by vector xp and the rotational axis of the planet. When the displacement L <C V^R, the 
effect of the translational quadrupole Mxp^Xp'' on the light deflection is relatively small, and the orientation of the 
quadrupolar pattern is defined by s as shown in subplots (e) in Figs. 6, 7. If i ~ \/J2R, the quadrupolar pattern may 
be complicated and hard to interpret as the deflection of light by the translational quadrupole Mx'ji^x-'p is comparable 
to that caused by J2- If, for any reason, ^/J-iR <C L, the deflection of light by the intrinsic quadrupole moment of 
the planet is swamped by the much stronger deflection caused by the translational quadrupole as shown in subplots 
(g) in Figs. 6, 7. 



C. Dynamic patterns 

Eqs. (51)-(53) can be used to work out a dynamic visualization of the deflection patterns as a function of time 
reckoned at the point of observation. Let us assume that a planet is traversing on the celestial sphere below a star 
along a great circle with a minimal separation angle Xm- As the deflection angles are rather small we can project 
this motion onto the plane of the sky so that it represents a straight line. We further assume for simplicity that this 
straight line is parallel and sufficiently close to the ecliptic. This is a good approximation since the orbital inclinations 
of Jupiter and Saturn are small and do not exceed 2.5° [82]. Eqs. (51)-(53) can be viewed as a one-to-one mapping 
of the planetary position specified by coordinates (X, Y) in the plane of the sky, to the locus of the star image shifted 
by the planetary gravitational field, at coordinates (x, y) (sec Fig. 8 and reference [83] for a more detailed geometry 
of the mapping). The undeflected position of the star (when the planet is at "infinity") is the center of the mapping 
which is also the origin of both coordinate frames: {X, Y) and {x, y) (notice that we have shifted Y axis for making 
better graphical visualization). In other words, the undeflected position of the star is at the point a; = 0, y = 
where we also have X = 0, y = 0. The planet moves along a straight line parallel to Y axis at a constant separation 
X = Xq from the undeflected position of the star. When the planet is at position angle <^ in the coordinates {X, Y) , 
the deflected position of the star is at angle 6 in the coordinates {x,y). Notice that in general, ^ ip because the 
dipolar and quadrupolar light-ray deflections are not purely radial. Precise correspondence between the two angles is 
established by the central mapping transformation equations. The angle ip runs from — 7r/2 to -|-7r/2 while the angle 
6 runs in the most common case from to 2-k (see Fig. 8 and the text below). 

Substituting the deflection angle cJm = (a;,y), sin — Y' I c?, cos ip — X/cl d = \/X^ + r2 to Eq. (51), and 
approximating for simplicity cosx = 1, we draw the monopolar deflection mapping equations in the following form 

- = 2r^^, y = 2r^^^, (54) 

where r = 2M/X, so that 2r is the maximal deflection angle reached at the time when the angular distance between 
the planet and the star is minimal (Xm)- The mapping given by Eq. (54) is purely radial which means that the 

deflected position of the star and the planet lie on a straight line passing through the center, and 9 = (p (see Fig. 8). 
As the planet moves from ip = —tt/2 to (p = +Tr/2 along the Y axis, the deflected position of the star describes the 
circle 

{x - rf +y'^=r'^ , (55) 

with radius r and the origin located at the point {x — r,y = 0). In polar coordinates, x = pcosO, y = psmO, and the 
correspondence between 6 and (fi is 6 = ip, as we have pointed out earlier. The monopolar deflection curve in polar 
coordinates is 

p = 2r cos ip , (56) 

which yields the circle shown in the left graph of Fig. 9, where x = A/3, y = AX. 

The dynamic curve of the dipolar deflection can be easier represented in coordinates {x',y') that are rotated 
clockwise with respect to the {x, y) coordinates through angle ip, where tp is deflned by Eq. (49, 

x' ~ x cos — y sin ip , y' = x sin ip + y cos ip . (57) 

In the new coordinates, the dipolar mapping transformation is deduced from Eq. (52) and is given by 

, x^(x^ — y^) x^Y 
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where p = rL/X. The dipolar mapping generally contains an orthoradial deflection component as demonstrated in 
Fig. 8. As the planet moves from f = —tt/2 to (f = +n/2 along Y axis, the deflected position of the star describes 
the curve 

shown in the middle graph of Fig. 9. This curve, called cardioid [83], is symmetric around the axis x' with a cusp at 

the origin as shown in the middle graph of Fig. 9. In polar coordinates, x' = pcosO', y = psm6', the correspondence 
between 0' and ip is 0' = 2ip, and in these coordinates the cardioid curve takes the form 

p = p{l + cos2ip) . (60) 

Its orientation depends on the angle tjj, and the magnitude on the ratio L/X. Notice that the vector xp was assumed 
to be constant. In reality, the planet may have a residual transverse velocity with respect to its ephemeris trajectory. 
In this case, the dipolar deflection curve will slightly deviate from the cardioid. 

The dynamic trajectory of the quadrupolar deflection becomes simple in coordinates (i, y) rotated clockwise with 
respect to the (a;, y) coordinates through angle 2tp, 

X = xcos2ijj — ysm2tjj , y = xsm2ip + y cos2il) . (61) 

In these new coordinates, the quadrupolar mapping transformation may now be inferred from Eq. (53) 

X^iX^-SV^) „ X^Y(SX^-Y^) , , 

^ = 4? (x2 + y2)3 ' y = 49 (x2 + y2)3 ' (62) 

where the parameter q = = {r/2){L/X)'^ for the deflection component caused by the translational quadrupole 

Mxp'xp'*, and q = q.j = {r/2)J2{R/X)^ for the term generated by the intrinsic quadrupolar oblateness J2 of the 
planet. The quadrupolar mapping contains an orthoradial deflection besides the radial component. As the planet 
moves from ip = —■jt/2 to (p = +7r/2 along Y axis, the deflected position of the star describes the curve 

4{x^ +f- qxf = 27q^{x'' + ff , (63) 

which is shown in the right graph of Fig. 9. This curve, called the Caley's sextic [83], is symmetric around the axis x 
and has a small secondary loop near the origin. The polar equation of the sextic is gained by substitution of the polar 
coordinates x — pcosO, y = psinO in Eq. 63, and solving the cubic equation with respect to p. The quadrupolar 
mapping correspondence between the angles 6 and (fi is 6 = 3(p, and the deflection caused by the translational 
quadrupole 

p = ?!, (cos3<^ + 3cos<^) . (64) 

The sextic curve for the deflection caused by the oblatness J2, takes the same form after replacing q^ qj and tp ^ u). 
The quadrupolar deflection curve has a secondary loop passing through the origin of coordinate axes as shown in Fig. 
9. The cross-over point of the secondary loop is 1/8 of the maximum deflection p = 4g'i:, in the mid-point of the 
primary loop. This means that the magnitude of the secondary loop is about few microarcseconds for Jupiter that is 
withing the range of experimental measurement by SIM and SKA. 

Fig. 10 shows the combined deflection trajectory from the dipole cardioid and the quadrupole sextic presented in Fig. 
9. The resulting curve looks similar to a cardioid of a somewhat larger amplitude, and is inclined at some intermediate 
angle with respect to the cardioid and sextic symmetry axes. Obviously, this curve swamps the quadrupolar deflection 
and can be confused with a pure dipolar one, because the vector xp is not known beforehand. This demonstrates 
the crucial importance of accurate determination of the spurious dipole component in any experiment designed to 
measure the intrinsic gravitational quadrupole ,/2. Since the ephemeris error is expected to be a smooth function 
of time, microarcsecond observations of multiple stars at different position angles to the direction of the dipole axis 
should help to separate these two contributions to the total deflection pattern. 



VII. EXPERIMENTAL TESTS OF RELATIVITY AND REFERENCE FRAMES 

Crosta and Mignard [30] assumed that as light propagates toward the observer the gravitating planet's center of 
mass is at rest at the origin of the chosen intermediate coordinate system in the sky. Our formalism stems from a 
more general dynamical interpretation, because it takes into account the effect of planet's motion on the deflection 
angle in a rigorous relativistic paradigm. The planetary motion appears in the light-ray deflection measurements in 
two ways: 
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• explicitly - in the form of velocity-dependent term, vp, and acceleration, ap, of the planet with respect to the 

coordinate system of our choice; 

• implicitly - in the form of the retarded time, Si = ti —ri, in the coordinates, xp = xp{si) of the planet as well 
as in its velocity and acceleration. 

General relativity may hypothetically be violated in the higher-order terms of approximation, which will become 
detectable in the near future by astrometric facilities operating at the 1 /zas level of precision. This presumable 
violation of general relativity can be parametrized and tested in light deflection experiments. Verification of general 
relativity and the quest for beyond-GR phenomena include the search for scalar fields, observation of gravimagnetic 
and retardation-of-gravity effects associated with the motion of planet, detection and measurement of the dipolar and 
quadrupolar anisotropy, and improvement of fundamental reference frames and planetary ephemerides in the Solar 
system. In what follows, we discuss these experimental opportunites in more detail by analysing different modes of 
the gravitational light-ray deflection. 

A. The monopolar deflection 

The monopolar deflection of light is given by Eqs. (42), (45), (46). It is standard experimental practice to introduce 
the PPN papameter 7 leading to the replacement [25] 

a ^ = — - — a (65) 

in Eq. (42) as well as Eqs. (43), (44). The parameter 7 characterizes deviation of gravity from a pure geometry 
and is associated with the presence of hypothetical scalar fields remaining from the epoch of the Big Bang [84]. This 
parameter is usually measured in the solar light-ray deflection experiments with a precision approaching to the level 

of 10^''^ 10-^ [4, 85]. 

Orbital motion of a planet generates hoi component of the gravimagnetic field which depends on the time derivative 
of the dipole moment as appears in Eq. (13). This gravimagnetic component leads to the explicit dependence of the 
monopolar light-ray deflection on the radial component of the orbital velocity ilp of the planet, effectively modulating 
the planetary mass as shown in Eq. (46). The gravimagnetic term was first derived in [9] and analytically confirmed 
in [69, 86, 87]. It can reach about 1.6 /xas for a star observed at the limb of the Jupiter. 

The gravimagnetic term is the result of the asymmetric exchange of energy between the photon and the gravitational 
field of the moving body. It is well-known [48] that in case of a static gravitational field the overall angle of light 
deflection does not depend on the energy (frequency) of the photon. The photon gains energy as it approaches the 
gravitating body, and loses it as it moves away. If the body is at rest with respect to the chosen coordinate system, 
the gain and the loss of energy cancel out, but this balance is violated if the body is moving radially. Appearance of 
the Doppler modulation of the deflection of light can be also understood if one remembers that the deflection angle 
is associated with spatial components of the null vector of photon. Lorentz transformation from a static frame of the 
planet to the moving frame of observer transforms and mixes all four components of the wave vector and its time 
component affects the spatial components in the moving frame. This is clearly shown in the paper by Klioner [69]. 

The gravimagnetic deflection of light can be measured in a single-epoch observation only if (1) mass of the planet is 
known sufficiently well; (2) the parameter 7 is excluded since it correlates linearly with the gravimagnetic term. On the 
other hand, if we can conduct observations in different epoches, the gravimagnetic term in Eq. (46) will periodically 
modualte the radial deflection because of the relative motion of the planet with respect to observer. Hence, it can be 
measured independently of the parameter 7. In case of the Jupiter the amplitude of this gravimagnetic modulation 
reaches 3.2 /xas and has a main period of one year due to the orbital motion of the Earth. 

B. The dipolar deflection 

The dipolar deflection of light is given by Eqs. (43), (45) where the coordinates of the planetary center of mass 
Xp = xp{si) depend on the retarded instant of time si = ti — ri, which is a solution of the retarded null-cone 
equation (15). These retarded coordinates of the planet deflne the dipole moment /(.si) = A/.?p(.si) in the multipolar 
expansion of the gravitational fleld of the planet. In general relativity, the law of conservation of linear momentum 
allows us to eliminate the dipole moment by placing the center of mass of the planet at the origin of the coordinate 
system. 

If xp{si) = 0, all coordinate-dependent effects in the gravitational deflection of light vanish. It means that the 
planet deflects light from the retarded position defined by the retarded solution of the gravity wave equation (11) 
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where t = ti is the time of observation and r = ri is the distance between the observer and the planet as shown in 
Fig. 3. This retardation is due to the fact that hght and gravity interact on the hypersurfacc which is the intersection 
of two nuU cones - the past hght cone of observer and the future gravity cone of the moving planet. This important 
prediction of general relativity [20] was tested and confirmed in 2002 [21-23] with precision of ~ 20%. SIM, Gaia, and 
SKA can improve this measurement by, at least, a factor of 10. In a practical experimental setup, the retardation of 
gravity effect can be parametrized by introducing a parameter e to the retarded time, s = t—er [88] in the dipolar 
light-ray deflection angle (43). If e = 0, the gravitational field of the planet would deflect photons instantaneosly. 
The general relativity prediction is that e = 1, which corresponds to the case of gravity propagating with the same 
speed as light. This consideration and comparision with the calculations of Crosta and Mignard [30] make it evident 
that the position of the Jupiter's center of mass, that is not specified in the paper [30], must be taken at the retarded 
instant of time Si = ti — ri with respect to observer. 

For arbitrary e 7^ (that is when the speed of gravity and light differ) light is deflected by the planet from its 
orbital position taken at the retarded time s^. It leads to a non- vanishing displacement vector L 

L = xp{si) - xp{se) = (e - l)vpri - ^(e - ifaprj + 0[{e- if] , (66) 

where vp and ap are velocity and acceleration of the planet respectively. If gravity operates on a null cone that is 
different from the light cone, it brings about the dipolar anisotropy 

do = (e — l)a — {vp ■ n)n — {vp ■ rn)m (67) 

— -{e—lfa-j {ap ■ n)n^ — {ap ■ in)m^ 

which can be measured to evaluate the retardation of gravity parameter e [22, 23, 26, 90]. For a star observed 

on the limb of Jupiter {d ^ R where R is the radius of the planet) the velocity-dependent term in Eq. (67) is 
^ <^iimh{v / c){r I / R) ~ 8140 ytias, and the acceleration-dependent term is ~ aiinih{v^ /c^){ri/2R) ~ 0.2 y^as. The reader 
should notice that the dipolar deflection (67) depends on the transversal component of the planetary velocity only. 
Since the velocity of the planet is almost lying in the ecliptic plane, it makes the position angle tp = — 7r/2 of the 
unit vector z = L/L in Eq. (49) so that the symmetry axis of the cardioid in Fig. 9 gets aligned with Y axis. This 
also allows us to physically interpret the dipolar deflection associated with the retardation of gravity as caused by the 
gravimagnetic field of the Jupiter originating from its translational orbital motion [27-29, 89, 90]. 

A VLBI experiment was conducted in 2002 [21] to search for the residual dipolar deflection of light from the QSO 
J0842-I-1835 by the moving gravitational field of the Jupiter, caused by the presumable difference between the speed 
of gravity and light parametrized by e. This experiment did not show any deviation from the general relativistic 
model of light deflection by the moving planet in excess of 20%. 

Assuming that general relativity is fully valid, measurements of the dipolar anisotropy of light deflection can be 
used to determine the position of the center of mass of the planet on its orbit more accurately than all other currently 
available astrometric techniques (see section VIII A for numerical estimates). 



C. The quadrupoleir deflection 

It is evident from Eq. (44) that in order to measure the quadrupolar deflection of light caused by the planetary 
oblatness J2, the following condition on the dispacement vector xp must be satisfied 

\xp\ < \aL\ < ^/J^R , (68) 

where ctl is the measurement error of xp. If condition (68) is fulfilled, all terms in the second line of Eq. (44) 
depending on the ratio {L/df can be abolished. In fact, cri, 7^ because the planetary ephemeris is known with 
a limited accuracy. Hence, one should consider the vector xp as a fitting parameter that is to be determined from 
light-ray deflection observation by successive iterations. JPL ephemerides provides the initial value of xp for such 
iterations, which can be improved if the measurement accuracy on the deflected positions is high enough (see next 
section). For the Jupiter \/J2R — 8500fcm, so that in order to measure the intrinsic quadrupole moment J2 of the 
planet its orbital position must be known to < 8500fcm. The currently available JPL ephemerides for major Solar 
system bodies, determined from direct optical/radio observations and spacecraft tracking, are believed to be more 
accurate. However, astrometry of gravitational light deflection provides another independent method of reflning the 
dynamic model of the Solar system. A non-zero value of L also emerges if the time of observation ti instead of the 
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retarded time si = ti — ri is used in the light-ray deflection equations whenever the instantaneous position of the 
planet is computed, which would be the case of e = in Eq. (66). Then, L w 35000fcm for the Jupiter, and condition 
(68) is violated, thus, making the intrinsic quadrupolar deflection of light unmeasurable. 

After the displacement xp is determined from the estimated dipole anisotropy, one can shift the origin of the 
coordinate system in the sky to the center-of-mass of the light-ray deflecting planet to supress the spurious dipolar 
term as much as possible. If general relativity is a correct theory of gravity, this will reduce both the dipolar and 
quadrupolar anisotropies caused by the translational displacement L of the planetary center of mass. On the other 
hand, if general relativity is violated one may not be able to remove the dipole anisotropy because the structure of the 
light-ray deflection equations will be different from those shown in Eqs. (42)-(44). This would make it more difficult 
to measure J2 in the models of the light deflection based on an alternative theory of gravity. 

VIII. OBSERVATIONAL CAPABILITIES 

Three astromctric missions capable of angular measurements at about 1 /las level accuracy are under development 
for the relatively near future. Two of them, the Space Interferometry Mission (SIM) [11] and Gaia [14], are space- 
based optical facilities, while the Square Kilometer Array (SKA) [12] is a ground-based radio telescope array. We 
discuss in more detail the SIM and SKA since the gravitational applications of Gaia has been described in [30]. 

A. SIM optical interferometry 

Perhaps, the SIM, which is a Michelson-type interferometer with articulating siderostat mirrors, holds the best 
prospects for precision tests of general relativity in the Solar system through gravitational bending effects. In these 
experimental-gravity applications the adavantages of the SIM facility are as follows: 

1. unlike Gaia, SIM is a pointing mission which can be aimed at an interesting object at the most appropriate 
time. Planc;tary near-limb grazing passages of sufficiently bright stars will be necessary to confidently measure 
the monopolc light-ray deflections from all the major planets, and the dipole and quadrupole moments from the 
Jupiter and, presumably, the Saturn. Such conjunction events are fairly rare but can be accurately predicted. 

2. SIM is primarily designed for astrometric detection of terrestrial-mass planets. In the differential regime of 
operation the interferometer is expected to achieve the unprecedented accuracy of 1 /xas in a single observation 
on stars separated in the sky within ^^2 degrees. A single measurement accuracy of Gaia is more than an order 
of magnitude poorer even for the brightest stars. 

3. the baseline of SIM can be rotated through 90 degrees for a dedicated observation, providing nearly simul- 
taneous two-dimensional observations of a set of stars while the light-ray deflecting planet is still near them. 
Two-dimentional observations on a given set of stars are crucial for unambigious disentangling the dipole and 
quadrupole deflection patterns, which may have similar magnitudes in the vicinity of a planetary limb. 

4. SIM will self-calibrate its 15°-wide field of regard to rms 20-40 pm in Zernike polynomial coefficients, which 
corresponds to the range of 0.4 - 0.9 /las subtended on the sky. This dramatically reduces the correlated and/or 
systematic errors and allows us to extend precision of differential measurements at 1 /uas to wider angles in the 
central part of the field of regard, if needed. 

5. owing to the small field stop in the optical design, SIM can observe stars and quasars as close as several 
arcseconds from the planetary limb, that is for impact parameter of the light ray comparable with the planetary 
radius. 

Assuming that four stars are observed with SIM at roughly the same impact angular distance x from Jupiter's center 
of mass, the signal-to-noise ratio (S/N) in measurement of the dipole deflection angle ud from a single one-dimensional 
observational visit is 

'S_ 

where ctq is the single measurement differential error, and (/ < 2) is a geometric factor that depends on the actual 
position of the stars and the angle between the baseline of the SIM interferometer and the direction of vector xp 



0.0254/ 



1 /ias 




1' 




z 






.X. 




1 km 



(69) 



19 



characterizing the dipolar hght-deflection pattern (see Fig. 1). The measurement error gl on the magnitude of 
L = \xp\ for this configuration is 

2 

km. (70) 

Thus, the best achievable accuracy in determination of L (the error in the position of the planetary center of mass) 
with SIM is roughly 20 km in the sky plane, when x = 1' and co = 1 /ias. This is more than 10 times better than the 
currently available positional ephcmcris for the Jupiter [91-93]. However, it will be hardly possible to achieve this 
accuracy in a single measurement, because of the difficulty to find simultaneously four stars within 1' of the planetary 
limb. Multiple observations of diff'erent star configurations should be planned. 

Assuming that the radial distance between the observer and the Jupiter is 5 AU, the amount of the gravitational 
quadrupole deflection of light from the planet is roughly 

3 

/las , (71) 

as follows from the J2-dependent term in Eq. (44). It will be accurately measured with SIM if fairly dense fields 
of background stars can be found along the path of the planet in the sky. For the previously considered idealized 
configuration of four equally-distant stars, two one-dimesional observations with orthogonal baselines will be required 
to disentangle the dipolar and quadrupolar components of the light-ray deflection pattern. For a pair of mutually 
orthogonal observations, the signal-to- noise ratio in measurement of the quadrupole deflection angle ag is 
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where again (/ < 2) is a geometric factor. Interestingly, the quadrupolar deflection from the Jupiter is comparable 
to the expected dipolar deflection at x = 1'; but since ag falls off faster than an as the angular distance x of a, 
star from the planet increases, multiple observations at wider angles from the planetray limb may be advantageous to 
separate and to determine the dipolar component of the light-ray deflection. Because of the complicated pattern of 
the quadrupolar component, a combined reduction of all observed events will be required, after the dipolar correction 
will have been determined and taken into account. 



B. SKA radio interferometry 



Extensive discussion of various fascinating science drivers and of the evolving technical possibilities has led to a 
concept for the Square Kilometre Array (SKA) and a set of design goals [12]. The SKA will be an interferometric 
array of individual antenna stations, synthesizing an aperture with a diameter of up to several 1000 kilometers. A 
number of configurations will distribute the 1 million square meters of collecting area. These include 150 stations each 
with the collecting area of a 90 m telescope and 30 stations each with the collecting area equivalent to a 200 meters 
diameter telescope. The sensitivity and versatility of SKA can provide ~ 1 /zas astrometric precision and high quality 
milliarcscc-resolution images by simultaneously detecting calibrator sources near the target source if an appreciable 
component of SKA is contained in elements which are more than 1000 km from the core SKA [13]. 

Measurement of the light bending by a moving planet with microarcsecond accuracy requires a continuous phase- 
referencing observation of the target and the calibrating radio sources [21, 94]. The main limitation of the accuracy is 
the tropospheric refraction which affects radio observations. The large-scale tropospheric refraction can be estimated 
by observing many radio sources over the sky in a short period of time. At present the determination of the global 
troposphere properties can only be estimated in about one hour, and smaller angular-scale variations can not be 
determined in most cases. However, the SKA, by using observations in ten subarrays, on strong radio sources around 
the sky, will determine the tropospheric properties on time-scales which may be as short at five minutes. 

Quasars as astrometric calibrators have one peculiar property: they are variable. The radio emission from the 
quasar comes from the base of the jet which is formed by an accretion disk around a massive black hole. The massive 
outflows and shocks in the jet change the intensity and the structure of the radio emission. Hence, the position of the 
quasar reference point (a radio core) is variable by about 0.05 mas in most quasars [95]. Thus, the calibrators used 
to determine the SKA astrometric precision to better than 10 jisec have jitter which is somewhat larger. In order to 
reach the intended angular precision, the change in position of the calibrators must be determined. With knowledge 
of the evolution of radio cores and the monitoring of a basic set of primary calibrators the calibrator grid can be 
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determined, at least to 10 iias or better [13]. The use of many calibrators in the field of view would also diminish the 
net effect of position jitter. 

In addition to various special and general relativistic effects in the time of propagation of electromagnetic waves 
from the quasar to the SKA-VLBI antenna network, we must account for the effects produced by the planetary 
magnetosphere [96, 97]. To be more specific, let us concentrate on the case of the Jupiter. Consideration of the 
magnetospheric effects for the Saturn can be done similarly. 

Measurements obtained during the occultations of Galileo by Jupiter indicate [98] that near the surface of Jupiter 
the electron plasma density reaches the peak intensity A^o = 1-0 x 10^*^ m~^. We shall assume that the jovian magne- 
tosphere is spherical, though in reality the magnetosphere has a dipole structure and our model may underestimate 
the plasma content along the polar directions. We postulate that a radial drop-off of the plasma density N{r) is 
proportional to l/r^^"^ where r is the distance from the center of Jupiter, and A is a number. The guess is that 
A > 0, and we will assume that A ^ for the worst possible case. Hence, radial dependence of the electron plasma 
density is taken as N{r) = iVo(i?j/r)2+^, where Rj = 7.1 x 10^ m is the mean radius of Jupiter. 

The plasma produces a delay AT in the time of propagation of radio signal which is proportional to the column 
plasma density in the line of sight given by integral [99] 

^ ^ r" N{r)dr n N{r)dr 

where and ri are radial distances of quasar and radio antenna from Jupiter respectively, and d = j^j is the impact 
parameter of the light ray from the quasar to Jupiter. In the experiment under discussion the impact parameter is 
much less than both ro and ri . Hence, 

^,(m-) = ^.ii.(^)"""^^, (74) 

where r(2;) is the Euler gamma-function. The plasma time delay (in seconds) 

AT(s) = AOAc-^v-'^Ni , (75) 

where c is the speed of light in vacuum measured in m/sec, v is the frequency of electromagnetic signal measured in 
Hz. 

It is worthwhile noting that, in fact, the SKA-VLBI array measures difference in path length between the radio 
telescopes. Hence, one has to differentiate A^"; in expression (75) with respect to the impact parameter d and project 
the result on the plane of the sky. This gives a magnetospheric VLBI time delay of 

Amag(ps) = 6.3 X 1Q-\A + 1)^ (^)' ^ . (76) 

normalized to the frequency = 8.0 GHz. The deflection of light Sumag. caused by the magnetosphere is evaluated 

as Samag. — cA,„ag/& rad. 

For example, substituting d = 5Rj and taking the SKA-VLBI baseline b — 6000 km we find 

<5Q!r„ag. - 163(z^o/j^)%as (^ = 0), (77) 

fon,ag. ^ 24 {i^o/i^f A^as (A = 1) , (78) 

(5a„ag. ^ 14 {lyo/iyf /ias (^ = 2) . (79) 



This represents the pure bending from Jupiter's magnetosphere at a single radio frequency Vq = 8 GHz, which should 
be compared with the magnitude of the gravitational deflection of light by various multipoles at the light-ray closest 
approach. SKA higher frequency is expected to be 35 GHz [100]. At this frequency the magnetospheric deflection is 
smaller as compared to 8 GHz by a factor of 20. 

In any case, the magnetospheric deflection estimate reveals that a single frequency observation of the light deflection 
will be affected by the magnetosphere at the level exceeding 1 ^as. This assumes that we should observe at two 
widely spaced frequencies to determine and eliminate the magnetospheric effects. However, to keep sensitivity at 
higher frequencies two polarizations and rather wide bandwidth must be used. The noise due to turbulence in the 
magnetosphere (and the Earth ionosphere) may also be a limit. However, this rapidly fluctuation model is fairly 
pessimistic and unlikely, and would probably average out to the steady state model. The most optimum method to 
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deal with the possible jovian magnetosphere component is not yet known and should be a matter of special study in 
radio astronomical community. 
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FIG. 1: Light propagates from star to observer in the direction shown by the unit vector k. The planet is displaced from the 
origin of the coordinate system by vector L that is a parameter of the data analysis algorithm. Its rotation axis is specified 
by the unit vector s. The impact parameter of the light ray is d. Two unit vectors, n and m, are orthogonal to the vector k 
and form a "plane of the sky" that is perpendicular to the line of sight, x denotes the angle between the star and the assumed 
position of the planet in the coordinate origin. 
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FIG. 2: Gravitational field of a planet is a retarded solution of the gravity-field wave equation (ff). In general relativity, 
gravitational field progresses on the hypersurface of a null cone from the past to the future. Directions of the propagation of 
the gravitational field of the planet are shown by arrows. The picture assumes that the coordinate-dependent effects associated 
with the dipole moment P are excluded {P = 0). Observer measures the gravitational field at time t when the planet is located 
at the retarded position on its orbit at the retarded time s = t — r, which is a null characteristic of the gravitational field. 
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FIG. 3: Gravitational field of a planet (the case of dipole 7' = is shown) affects only the particles lying on the hypersurface 
of the future gravity null cone. A photon emited by a star at time to arrives to observer at time ti along a null direction of 
the past light cone with a vertex at the observer. Therefore, the future gravity null cone of the planet and the past light cone 
of the observer must coincide along the null direction that is a null characteristic of the retarded solution of the gravity-field 
wave equation (11). The photon detected at time ti, is defiected by planet's gravity force from the planet's retarded position 
taken at time si = ti — ri. This effect of the retardation of gravity can be observed by measuring the amount of gravitational 
deflection of light by a moving planet, and used to determine the speed of gravity field propagating from the planet to the 
point of observation (see section V and papers [21, 22]). Note that the retardation of gravity would not be measurable if the 
planet were at rest with respect to the observer. 
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FIG. 4; Light-ray deflection by a static (a) and moving (b) planet. In case (a) the distance between the planet and observer 
does not change as light propagates. Thus, measuring the deflection of light does not allow us to determine experimentally 
whether the gravity force of the planet acts with retardation from position (1), or instantaneously from position (2). In case 
(b) the distance between the planet and the observer varies as the photon travels toward observer. The retarded interaction of 
gravity with light becomes apparent since measuring the angle of the gravitational deflection of light allows us to distinguish 
between positions (1) and (2) of the planet on its world- line making the causal structure of the gravity cone clearly visible and 
measurable. 
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FIG. 5: Local tangential coordinate system in the plane of the sky related to the global ecliptic coordinate system is shown. The 
unit vectors e\ and e}? point in the direction of increasing ecliptic longitude and latitude, respectively. The angular distance 
of the undeflected position of the star from the origin of the coordinate system is x (cf- Fig. 1), which is at position angle ip 
from the north direction. The total deflection a has both radial, a„, and orthoradial, a™, components: d — anfi + amrn. 
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FIG. 6: Gravitational deflection of light by the multipolar fields of the Jupiter. The magnitude of deflection is coded by 
colors in a logarithmic scale shown at the top. The vectors in each plot indicate only the direction of displacement, not the 
magnitude. The monopole deflection pattern is shown in plot (a). The artifactual dipolar deflection patterns caused by a 
positional displacement xp = Lz of the planet's center of mass from the origin of the coordinate system used for multipolar 
expansion of the gravitational field, are depicted in the central row of plots, for Li = 350 (b), L2 — 3500 (c), and L3 = 35000 km 
(d). The bottom row of plots (e-g) shows the combined quadrupolar deflection patterns generated by the intrinsic oblateness 
of the planet and by the coordinate-dependent translational moment, for the same values of L. The direction of rotation axis s 
is indicated with a thin black line. Note that the quarupolar deflection at Li = 350 km is dominated by that from the intrinsic 
moment, whereas for L3 = 35000 km, on the contrary, the translational quadrupole exceeds the light deflection caused by J2. 
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FIG. 7: Deflection of light by the Saturn caused by its mass (a), dipole (middle raw), and quadrupole (bottom raw) moments. 
Magnitudes of the deflection are expressed by colors in the logarithmic color grade shown in the upper right corner, and the 
direction (not the magnitude) of the deflection is shown by small arrows. Vector displacement L = Lz of the center of mass of 
Saturn from the origin of the coordinate system takes values (from left to right) of Li — 250 (b), L2 = 2500 (c), and L3 = 25000 
km (d) . The unit vectors z and s indicate the directions of the displacement L and rotational axis of the planet respectively. 
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FIG. 8: Geometry of the central mapping used for computation of the dynamic patterns (curves) of the gravitational deflection 
of the apparent position of a star by a planet as a function of time. The origins of the two coordinate charts {x,y) and X,Y 
coincide and placed at the undeflected position of the star. For the sake of graphical convenience, Y axis is displaced along 
the positive direction of X axis at the constant distance Xo- The available planetary ephemeris predicts that planet moves 
continuously from left to right along Y axis with X fixed as time passes, and two positions of the planet are shown to explain 
the mapping geometry. The dashed line is parallel to Y axis and shows the trajectory of the real planetary center of mass 
which may differ from the planetary ephemeris based on the other (independent) set of observations giving rise to the dipolar 
component of the light deflection. The minimal distance between the planet and the undeflected position of the star is X = Xo- 
Apparent position of the star is deflected from its unperturbed (catalogue) position in both radial and orthoradial directions 
that are a„ and Um respectively. The total deflection a — a„n + amfh. The mapping establishes a one-to-one correspondence 
F between the planet's coordinates {X,Y) and the star's deflected position {x,y), that is (_F : d ^ a, — > 9). 
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FIG. 9; Trajectories of apparent positions of a star deflected by the Jupiter. Left, middle and right graphs show respectively 
the monopolar (circle), dipolar (cardiod), and quadrupolar (the Caley's sextic) light-ray deflection curves. These curves are 
obtained by central mapping transformation explaind in Fig. 8 and section VI C. The planet passes directly below the true 
position of the star at the coordinate origin, at a constant velocity. The impact parameter Xo = 40 arcsecond, and the orbital 
error is L = 250 km. The time step between individual points is 0.02 days for the monopole deflection curve, and 0.002 days 
for the other two. Dynamic light-ray deflection curves for the Saturn have smaller amplitude but similar shape. 



33 



Dipole + Quadrupole 



T 




A X \ias 

FIG. 10: Combined dipolar and quadrupolar deflection of light by the Jupiter for the same parameters and configuration as in 
Fig. 9. Details of the quadrupolar deflection (the sextic curve) are hidden by the dipolar deflection (the cardioid). Separation 
of the two deflection terms is a non-trivial experimental problem, and generally will require multiple observations of a number 
of sources around the planet. 



